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Abstract 

Background: In the 1970s, mycoviruses were identified that infected the edible mushroom Lentinula edodes 
(shiitake), but they were not regarded as causal agents for mushroom diseases. None of their genes has been 
sequenced. In this study, the dsRNA genome of a mycovirus recently found in a shiitake commercial strain was 
sequenced and its molecular structure was characterized. 

Methods: A cDNA library was constructed from a dsRNA purified from the fruiting body of L edodes. The virus was 
tentatively named L edodes mycovirus HKB (LeV). Based on the deduced RNA-dependent RNA polymerase (RdRp) 
sequence, phylogenetic analysis of LeV was conducted. Because no virion particles associated with the dsRNA were 
observed by electron microscopic observation, atomic force microscopy (AFM) observation was chosen for 
achieving molecular imaging of the virus. 

Results: The 1 1,282-bp genome of LeV was obtained. The genome encoded two open reading frames (ORFs). 
0RF1 coded for a hypothetical protein and 0RF2 for a putative RdRp, respectively. In addition, a region coding for 
a NUDIX domain was present in 0RF1. There was a 62-bp intergenic region between 0RF1 and RdRp. Similarity 
with coat protein of mycoviruses was not found within the whole sequence. Based on phylogenetic analysis of the 
putative RdRp sequence, LeV grouped into a clade with dsRNA found in the basidiomycetes Phlebiopsis gigonteo 
and Helicobosidium mompo. The clade was placed apart from the Totiviridae and Chrysoviridoe families. As 
suggested from the genome sequence, AFM revealed that the structure of LeV was linear unencapsidated dsRNA. 

Conclusions: The results suggest that LeV represents a novel family of mycoviruses, found thus far only among 
the basidiomycetes. 
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Background 

In the 1970s, viruses that infect the cultivated mushroom 
Lentinula edodes, or shiitake, were extensively studied in 
Japan [1-3], and three morphologically distinct viruses 
were detected by electron microscopy [1,3]. However, 
unlike La France disease of the white button mushroom 
[4,5], mycoviruses have not been associated with shiitake 
diseases because these mycoviruses have commonly been 
found in healthy fruiting bodies [1,3]. In the USA, 
dsRNAs have also been observed in shiitake strains, but 
these appeared to be latent [6]. In the 1970s, shiitake cul- 
tivation was performed by inoculating mycelium spawn 
on oak logs; however, this labor-intensive method was 
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gradually replaced by indoor cultivation using sawdust 
substrate supplemented with rice bran. Currently in 
Japan, about 75,016 t (82% of all shiitake) are produced 
indoors annually using bag cultures with a sawdust-based 
substrate [7]. In a bag culture, the shiitake mycelium is 
fully grown in the substrate until brown pigment is pro- 
duced outside and the substrate becomes stiff. Complete 
browning of the exterior surface of the substrate is an 
important marker that normal fruiting bodies will 
develop in the following stage of cultivation. However, 
abnormal symptoms are occasionally observed in bag cul- 
tures, such as the growth of white or fluffy mycelia on the 
surface of substrate, inadequate or imperfect substrate 
browning [8], and malformations of the fruiting body. 
These symptoms often result in serious economic losses. 
Whether or not these abnormalities are associated with 
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mycovirus is unknown. As a first step toward addressing 
that question, 46 shiitake isolates belonging to 11 com- 
mercial strains were examined for mycovirus infections 
[8]. As a result, dsRNA was found in two isolates; one 
showed imperfect browning and the other was asympto- 
matic. Agarose gel analysis showed that the isolate with 
imperfect browning contained several dsRNAs, but the 
asymptomatic isolate contained only a single dsRNA. In 
this study, the single dsRNA was tentatively designated 
as Lentinula edodes mycovirus HKB (LeV) and was 
sequenced. 

Methods 

dsRNA isolation 

Three fruiting bodies were disrupted in 60 ml of 0.1 M 
phosphate buffer, and the virus fraction was precipitated 
with 10% PEG 8000 and 0.15 M NaCl, as described pre- 
viously [9]. The PEG precipitate was suspended in TES 
(10 mM Tris-HCl, 1 mM EDTA, 0.15 M NaCl, pH 7.0), 
and total RNA was isolated using the QIAmp Viral RNA 
Mini Kit (Qiagen) according to the manufacturers instruc- 
tions. Then dsRNA was isolated from the viral RNA by 
DNase I (Promega) digestion for 30 min at 37°C, followed 
by SI nuclease (TaKaRa) digestion. The resulting dsRNA 
was concentrated in nuclease-free water by filtration 
through Ultrafree 0.5 100 K centrifugal filters (Millipore) 
several times to remove degraded nucleic acids and salts. 
Finally, dsRNA was purified with the RNeasy Mini Elute 
Cleanup Kit (Qiagen). 

cDNA library construction and sequencing of dsRNA 

The purified dsRNA served as a template for cDNA synth- 
esis by random priming with the PrimeScript 1 st strand 
cDNA Synthesis Kit (TaKaRa) according to the standard 
protocol, except that the denaturing condition was chan- 
ged from 65°C, 5 min to 98°C, 1 min. The resulting 
cDNAs were electrophoresed in an agarose gel, and 
cDNAs sized 1.5-2.0 kb were extracted. The cDNAs were 
blunt-ended with T4 DNA polymerase, ligated into 
pUC118/HincII/BAP (TaKaRa), and transformed into 
E. coli DH10B cells by electroporation (Gene Pulser, Bio- 
Rad). The resulting cloned DNA was sequenced with the 
BigDye Terminator v.3.1 Cycle Sequencing Kit (Applied 
Biosystems) with M13 forward and reverse primers. 
A contiguous region (contig) was assembled with 
Sequencher™ 4.6 (Gene Codes Corp). 

AFM microscopic observation 

The PEG precipitate was suspended in 500 [i\ of TE and 
filtered through a 0.2-|im filter (Millipore), and 1 \i\ of 
the filtrate was diluted with 50 [i\ TE + 10 mM MgCl 2 . 
The sample was observed using atomic force microscopy 
(AFM) as previously described [10]. A total of 10 [d of 



the RNA solution was dropped onto freshly cleaved mus- 
covite mica (lxl cm), which after standing for several 
minutes was washed with distilled water. The sample was 
dried under a stream of nitrogen. Observations were per- 
formed on a MFP-3D (Asylum Research) in the tapping 
mode in air, using a silicon cantilever OMCL-AC240TS 
(Olympus). Fields of 2 [im x 2 [im were scanned at a fre- 
quency of < 2 Hz. To confirm that the AFM image 
observed was RNA, AFM imaging was also performed 
with PEG precipitate after RNaseA digestion. The length 
of the RNA molecule was measured using software devel- 
oped by the Research Institute of Biomolecule Metrology 
Co., Ltd. (Japanese patent P2000-230823A). 

Phylogenetic analysis 

An unrooted polygenetic tree was constructed with 
sequences retrieved by a PSI-BLAST [11] search (Table 1) 
using the neighbor-joining method with a bootstrap of 
1000. The sequence identity was below 40%; thus, we used 
the Multiple Alignment using Fast Fourier Transform 
(MAFFT) program L9INS-I http://mafft.cbrc.jp/align- 
ment/server/[12] for making the multiple alignments and 
constructing the phylogenetic tree. Sequences with boot- 
strap values above 70% were included in the tree and 
visualized on the web http://www.genome.jp/tools/mafft/ 

Results and discussion 

Genetic analysis 

A cDNA library was constructed from dsRNA purified 
from the fruiting body of Lentinula edodes strain HKB. 
An 11,282-bp cDNA contig containing two open reading 
frames (ORFs) was obtained (deposited under Accession 
No.AB429556) (Figure 1). ORF1 encoded a 218,428-Da 
protein composed of 1,975 amino acids, containing the 
conserved motif of a NUDIX domain (between 319aa 
and 451aa) (pfam00293, NUDIX, 1.56e" 04 ) [21]. Viral 
coding NUDIX domain has been found only in Pox- 
viruses [22]. LeV, as well as Phlebiopsis gigantea mycov- 
irus dsRNA 1 is the first mycovirus described, that codes 
for NUDIX domain. ORF2 encoded a 162,240-Da protein 
composed of 1426 amino acids, containing the conserved 
motif of RNA-dependent RNA polymerase (RdRp) 
(pfam02123, RDRP_4> 4.48e" 22 ). The deduced amino acid 
sequence of ORF1 showed significant similarity (3e~ 83 ) to 
hypothetical protein PgV-l_gpl of Phlebiopsis gigantea 
mycovirus dsRNA 1 [13]. The deduced amino acid 
sequence of ORF2 showed similarity with members of 
Totivirus [23] and Chrysovirus [24]. The family Totiviri- 
dae includes all the viruses of fungi and protozoa 
that have virions and a single-component dsRNA gen- 
ome. The family Chrysovirus includes dsRNA viruses that 
infect fungi or plants and have four genome components 
in isometric virions. Because LeV contained a single 
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Table 1 RdRp sequences of viruses included in the phylogenetic tree 



Virus 


Accession no. 


E-value 


Identity 


Similarity 


References 


Lentinulo edodes mycovirus HKB LeVHKB 


DA 1111 (~\~7 








This study 


Phlebiopsis giganteo mycovirus dsRNAI PgV1 


t~ A I O A O O O 


/I Q -1 50 

4e 


~>ci~> /1 1 r\~7 l~)~\C\/\ 

383/1 1 97 (32%) 


580/1 1 97 (48%) 


[13] 


Helicobosidium mompo V670 L2-dsRNA virus HmV670 


a d ~>~7 r ~>oo 

Ad27528o 


7e 


26/65 (40%) 


40/65 (62%) 


unpublished 


Helminthosporium victoriae 145S virus HvV145S 


AAM6895i 


1 e 


89/328 (27%) 


1 a 1 /~>~>o //nn/\ 

1 41/328 (43%) 


unpublished 


Verticillium chrysogenum virus VcV 


a r^r~ inn 

AUG21213 


1 e 


99/406 (24%) 


i~7~>//ir\^" //nn/\ 

1 72/406 (42%) 


unpublished 


Cryphonectrio nitschkei chrysovirus 1 CnCV1 


A r — 1 — ?mr i 

AC 1 79257 


9e 


83/318 (26%) 


1 36/31 8 (43%) 


[14] 


rci IILIIIIUi 1 1 LiliybUytzilUill VllUb rLV 


A AK/IQ^m 
aaivi^jdu 1 


Dc 


i uy/^-j J \Z*-\- /0) 


1 / J/ T-JJ \DO /0) 




Anthurium mosaic-associated virus AmV 


ACU 11563 


Se-' 1 


67/227 (30%) 


101/227 (45%) 


unpublished 


Grapevine associated chrysovirus-1 GaCV1 


ADO60926 


5e -09 


86/371 (23%) 


138/371 (37%) 


[16] 


Aspergillus fumigotus chrysovirus AfCV 


CAX48749 


4e- 08 


102/455 (22%) 


1 74/455 (38%) 


[17] 


Circulifer tenellus virus 1 GTV1 


ADK12924 


5e- 07 


82/327 (25%) 


135/327 (41%) 


[18] 


Spissistilus festinus virus 1 SpFV1 


ADK12922 


1e -06 


84/338 (25%) 


139/338 (41%) 


[18] 


Soccharomyces cerevisioe virus L-A ScVLA 


AAA50508 


2e -os 


58/230 (25%) 


100/230 (43%) 


[19] 


Sphaeropsis sopinea RNA virus 1 SsRV1 


AAD11601 


1e-° 3 


48/202 (24%) 


81/202 (41%) 


[20] 



E-values, identity, and similarity values were obtained from PSI-BLAST searches of the deduced amino acid data of RdRp sequences 
Viruses that showed bootstrap values above 70% were included in the phylogenetic tree 



component dsRNA genome, it could be classified as a 
member of Totivirus [23]. But unlike Totivirus, there was 
no coding region for a coat protein within the LeV 
genome. 

AFM microscopic observation 

No virion or vesicle particles were detected by electron 
microscopy observation of the PEG precipitate of LeV. 
Purification of the virus by sucrose density-gradient 
ultracentrifugation was unsuccessful. When the PEG pre- 
cipitate was directly digested with RNase, the viral 
dsRNA could no longer be detected by agarose gel elec- 
trophoresis (data not shown). These experimental data 
suggested that the virus did not form a virion and there- 
fore was vulnerable to excess purification. Because the 
method is nondestructive, AFM was chosen for further 
observation [25]. The PEG precipitate was suspended in 
TE and directly observed under AFM after being filtered 
through a 0.2-mm filter (Millipore). Numerous linear 
particles were observed (Figure 2A). The structure (Fig- 
ure 2B) was very similar to the AFM images of dsRNA 
observed by Abels et al. and Vilfan et al. [26,27]. To con- 
firm that these structures were RNA, PEG precipitate 
after RNase digestion was also observed by AFM (Figure 



dsRNA: 11,282 bp 

NUDIX domain 

5 , 0RF1 1 » 975aa ORF2 (RdRp) 1,426 aa _ g , 

899 6856 6889 11169 

Figure 1 Diagrammatic representation of genome organization 
of LeV. The dsRNA of LeV encodes two ORFs. ORF1 is an unknown 
protein and has a region with high similarity with the NUDIX 
domain. ORF2 encodes the RNA-dependent RNA polymerase (RdRp). 

v / 



3B). After the RNase digestion, the linear materials were 
no longer present. Thus, LeV was revealed as a linear 
dsRNA. The length of the dsRNA was 3,539 nm (Figure 
3A), and its molecular weight was estimated (based on 1 
[im = ~3 kb dsRNA [27]) to be 11.8 kb. The estimated 
size was very close to that of the sequenced dsRNA 
(11,282 bp). 

Phylogenetic analysis 

In the current classification of viruses, unencapsidated 
mycoviruses are placed in either the Hypoviridae 
[28,29] or Endornaviridae [30,31] families. However, 
the putative RdRp sequence of LeV showed no similar- 
ity to the RdRp sequences of Hypovirus or Endorna- 
virus. An unrooted polygenetic tree was constructed 
with sequences retrieved by a PSI-BLAST [11] search 
(Table 1) using the neighbor-joining method with a 
bootstrap of 1000. Sequences with bootstrap values 
above 70% were included in the tree (Figure 4). This 
analysis grouped LeV with Phlebiopsis gigantea mycov- 
irus dsRNAI (E value; 4e" 150 ) (CAJ34333) and Helico- 
basidium mompa V670 L2-dsRNA virus (translated in 
this study from AB275288) (E value; 7e" 33 ) in a clade 
independent of Totiviridae and Chrysoviridae (Figure 
4). Both P. gigantea and H. mompa are basidiomycetes. 
No virion particles associated with P. gigantea mycov- 
irus dsRNA have been described [13]. Currently, 
whether H. mompa dsRNA forms virions or not is 
unknown. If H. mompa dsRNA was also unencapsi- 
dated, this would suggest that there is a novel family 
of mycoviruses so far found only in the basidiomycetes, 
having a monopartite dsRNA genome but do not pack- 
age into virions. 
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Figure 2 AFM images of a LeV. AFM images of PEG precipitate. (A) LeVs appear as numerous linear nucleic acids. (B) Image of a single LeV 
dsRNA. 



Hypothesis of unencapsidated virus in the 
basidiomycetes 

L. edodes is a white-rot basidiomycete and P. gigantea is 
also a basidiomycete that causes white rot in conifer 
logs and stumps [32,33]. Often, P. gigantea is isolated 
from bark beetle, as is Lentinula boryana, a fungus 
belonging to the same family as L. edodes [34]. Virion 
structure is necessary for a virus to infect and exit from 
the host cell and gain a greater probability of propaga- 
tion. Because bark beetles feed on basidiomycete fungi 
[34], the evolution of unencapsidated virus in white-rot 
basidiomycetes may have been achieved through their 
association with wood-feeding insects. Because their 
host can be efficiently transferred to new environments 
by bark beetles, the coat protein would no longer be 
necessary for virus propagation. Additionally, the con- 
served amino acids of putative RdRp molecules of unen- 
capsidated viruses isolated from plant-feeding insects, 



such as the alfalfa hopper and beet leafhopper [18], are 
significantly similar to LeV (Figure 4) (Table 1). This 
fact also supports the hypothesis that unencapsidated 
species of mycovirus present in the basidiomycetes 
might be closely associated with insects. 

Conclusions 

dsRNA found in a commercial strain of Lentinula 
edodes (shiitake) (designated as L. edodes mycovirus 
HKB; LeV) was sequenced. The 11.8-kb genome con- 
tained two ORFs. ORF1 coded for a hypothetical protein 
containing a NUDIX domain. Previously, viral NUDIX 
domain was found only in Poxviruses. LeV is the only 
mycovirus that codes NUDIX domain. ORF2 coded for 
RdRp with high similarity with Totivirus and Chryso- 
virus. The genome did not code for a coat protein, and 
AFM observation revealed LeV to be unencapsidated 
linear dsRNA. The results show that LeV is a novel 
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LeVHKB 



new family? 




Chrysoviridae AmV 




Unclassified 
dsRNA virus 



Figure 4 Phylogenetic tree generated from the deduced amino acid sequences of RdRp. The unrooted phylogenetic tree of the aligned 
sequences was constructed based on the neighbor-joining method using 1000 bootstrap replications. Numbers at the node indicate percentage 
of bootstrap supports. LeVHKB: Lentinula edodes mycovirus HKB; PgV: Phlebiopsis gigantea mycovirus dsRNAI; HmVL2: Helicobasidium mompa 
V670 L2-dsRNA virus (translated in this study from AB275288); GTV1: Circulifer tenellus virus 1; SpFV1: Spissistilus festinus virus 1; ScVLA: 
Saccharomyces cerevisiae virus L-A; SsRV1: Sphaeropsis sapinea RNA virus 1; AmV: Anthurium mosaic-associated virus; HvV145S: Helminthosporium 
victorioe 145S virus; CnCV1: Cryphonectrio nitschkei chrysovirus 1; VcV: Verticillium chrysogenum virus; PcV: Penicillium chrysogenum virus; AfCV: 
Aspergillus fumigotus chrysovirus; and GaCV1: Grapevine associated chrysovirus-1 . 



mycovirus with a monopartite dsRNA genome that is 
closely related to Totivirus, but that it does not form a 
virion particle and might represent a new class of 
mycovirus. 
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